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In this paper, the effects of interfacial tension between the sediment solid particle and liquid on 
the settlement of sediment floes are investigated. After a discussion of mechanical and physical 
chemistry, we give a settling velocity expression including such dynamical information of the floe 
growth as interfacial tension and primary particle size etc.. The resulting expression indicates the 
average settling velocity of sediment floes increases with increasing solid-liquid interfacial tension in 
a form of power law and deceases with the primary particle size. We report on a general method for 
analyzing settling behaviors of sediment floes under different flocculation conditions and verify the 
rationality of the assumption of tension-induced flocculation by fltting typical experimental data to 
the electrolyte concentration-dependent sedimentation model which can follow from the relationship 
between interfacial tension and electrolyte concentration. 
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I. INTRODUCTION 

The aggregation-sedimentation of sediment floes, 
which is widely involved by many areas of industry and 
engineering such as mineral processing, irrigation works, 
and hydraulic engineering, is a basic phenomenon in the 
nature. Many fundamental researches have been devoted 
to understanding of it 
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ticles aggregate and form clusters or microflocs which 
may continue to combine the remaining single particles 
or other microflocs and thus become larger floes. When 
the growth of floes is induced by the Brownian motion, 
the diffusion limited cluster aggregation (DLCA) model 
or the reaction limited cluster aggregation (RLCA) model 
can well describe the process of this aggregation. How- 
ever, when the floes become large enough to settle under 
gravity, sedimentation alters the growth mechanism. If 
the sediment concentration is low, the floes fall under 
gravity, grow, fracturing and restructuring; the size of the 
floes reach a final value and the space-filling network con- 
necting does not take place. Conversely, if is sufliciently 
large, the growing microflocs will form the space-filling 
networks which spans all the cell very quickly. 

The DLVO theory HH, which attributes interactions 
between colloidal suspensions to Van der Waals forces 
and electrostatic repulsion, dynamically provides a quan- 
titative description of many facts of flocculation. In the 
DLVO framework, colloidal suspensions tend to coagu- 
late due to Van der Waals attractive forces, while elec- 
trostatic repulsive forces promote stability of the sys- 
tem, and the repulsive barrier may be decreased through 
changing the physicochemical conditions, temperature, 
ionic strength, PH, solvent quality and so on, which in- 
duces the flocculation of colloidal particles. Although the 
DLVO theory can describe basic facts of flocculation of 
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colloids, it still can not completely elucidate the cause 
of flocculation phenomena. For example, the DLVO the- 
ory encounter many difficulties in aggregation of charged 
colloidal particles in aqueous media such as fluid-crystal 
coexistence phenomenon j and void formation 

0,01, etc, where a long-ranged attraction is required. 
In addition, the DLVO theory is also reluctant to explain 
the transition from RLCA to DLCA 14], since the fact 
that high electrolyte concentrations lead to a fractal di- 
mension arising from DLCA model indicates that Van 
der Waals force cannot act as such a role that may make 
formation of the floe attributed only to pure Brownian 
collisions which means the binding force between floes 
should be independent of the floe itself. 

To overcome discrepancies between theoretic predic- 
tions and experimental observations, different versions of 
modiflcation to the DLVO theory, some of which have 
additionally considered non-DLVO interactions resul ting 
from acid-base 0|, steric 0|, and hydrodynamic |l7j 
interfacial forces, have been presented. However, most 
investigations into the effects of the interface are con- 
centrated on interactions between particles 0, 0, |23| . 
For a macroscopic system of sediment floes, whose sizes 
can often reach 500'-~^1000 /im, effects of interfacial ten- 
sion between the solid particle and liquid on mechanical 
balance of the floe are more necessary to be taken into 
account, and on such a large scale the DLVO force be- 
tween particles becomes negligibly small in comparison 
with gravity and hydrodynamic forces acting on the floe, 
so that dynamical reason for sediment flocculation may 
be attributed mainly to surface tension on the floe. 

In this paper, we emphasize the effects of the inter- 
facial tension between solid particles and liquid on the 
settlement of sediment floes. Few works about settling 
dynamics of sediment floes are involved in the growth 
mechanism of floes, and hence the aim of the paper is 
to introduce such dynamical information as interfacial 
tension, primary particle size (pps) etc. in the settling 
velocity expression of sediment floes. We present a gen- 
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eral method for investigating sedimentation of sediment 
floes under different flocculation conditions and discuss 
settling behaviors influenced by electrolyte concentration 
in relation with typical experimental data. 



II. THE SETTLING VELOCITY 

We now consider a system of sediment particles of the 
identical size in still water, and focus on such the situ- 
ation that the sediment concentration <f> is low so that 
spatially networking does not happen. After a process of 
aggregating, large floes fall under gravity, then average 
settling velocity of the system is written as 



(1) 



where tpi denotes the weight of the i-th group of floes 
with the same size Di , and v{Di) is (constant) settling 
velocity of the floe belonging to the i-th group. It is 
necessarily emphasized that Di should be understood as 
diameter of sphere of the same volume as the floe includ- 
ing inside water. The isometric diameter is preferentially 
considered here since some of experiments of measuring 
settling velocity of sediment floes are based on weight 
change determination that directly connects with the iso- 
metric size. 

Stokes settling velocity formula can be employed, if 
the floe are not so large that effects of turbulence on 
it might be ignored, for expressing the size-dependent 
sedimentation velocity v{Di) of the floe [23|, namely 
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where pAi.L represents the density of the i-th floe and 
liquid, respectively, v is the viscosity of the liquid, g is 
acceleration of gravity, ps is density of the solid particle, 
and ti denotes the porosity of the i-th floe. Here the 
relation ''p^Zp'^ = 1 — is used. Generally, Eq. (2) is 
applicable for the case that the size of the settling floe 
does not exceed 150 pm. 

Investigations of floe structure [2^ I2M 01 have shown 
that the density (or porosity) of the floe may be evaluated 
from the following relationship: 



Dn 



(3) 



where ds denotes a modifying factor, Dc is the diameter 
of collision, which denotes the size of the smallest sphere 
that fully encompasses the floe, Dq is the diameter of the 
primary particle comprising the floe, and a is a constant. 
And the mass fractal dimension df can be scaled as df — 
3 — a. The number TV of primary particles within a floe 
can be expressed in terms of the collision diameter Dc 



and the isometric size £>, respectively, namely 

where Ai is the structure prefactor, and 
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as a result, the relation between Dc and D is given by 

(6) 

Due to Eq. (3) and Eq. (6), the sedimentation velocity 
v{Di) can be rewritten as 
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where A2 = d^ . The size-dependent velocity ex- 

pression Eq. (7) gives the information of mass fractal 
dimension d f , hence it is also be a theoretic equation for 
measuring d/ of settling floe ^2^. 



III. THE ADDITIONAL PRESSURE OF 
SURFACE OF THE FLOC 

Mechanism of growth of the floe affects sedimentation 
of it, hence we have to take into account effects of the 
mechanism of the growth on settling floe. It is pointed 
out in this paper that the interfacial tension 7 on the 
solid-liquid interface (or the Gibbs dividing surface) be- 
tween sediment solid particle and liquid plays a key role 
in aggregation and sedimentation of particles or floes. 
Obviously, the interfacial tension 7 balances with the ad- 
ditional pressure difference p between both sides of the 
interface. The additional pressure difference p, which is 
associated directly with stress which the floe is suffering, 
is an important quantity determining the structure and 
thus settlement of the floe. For a regular surface, p is in 
inverse proportion to radius of curvature of the surface, 
and then, for a floe surface with the fractal feature, what 
properties does it have? 

We assume, to expand the floe volume an increment 
dV , the additional pressure needs to overcome the surface 
work of the interface area adding dS , that is. 
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and then from Eq. (5) we have 
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dS = nPsDldN = {TTl3,/Do)d[{l-e)D^]; dV = -nD^dD, 

where Ps is the proportion factor denoting the ratio of 
solid surface contributing to external surface of the floe 
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to total solid surface, which is associated with the shape 
of the floe. Combining Eq. (3), (8), and (9), we derive 
the expression of the additional pressure difference, 
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It is clearly seen that the additional pressure of the floe 
surface is inversely proportional to the size of the floe 
since the fractal dimension df is generally less than 3, 
and it decays following the power-law form, which is a 
reflection of abnormal specific surface area of the floe, and 
also the additional pressure may connect with electrolyte 
concentration through 7 according to Gibbs adsorption 
law. 

The reaction force of the additional pressure, acting on 
the floe, provides a binding force that can stick particles 
together. However, the inverse dependence on the floe 
size of the binding force characterized by p implies that 
not all positions on the surface of a certain floe can com- 
bine the other specific floe due to the impact of external 
forces such as gravity, shear force, etc. on the specific 
floe. In other words, under the given mechanical condi- 
tions, a floe would have a eombinable region relative to 
the other floe colliding with it, whcrc^ both of them may 
satisfy the mechanical condition of balance, as a result, 
a flrm agglomeration has to require many collisions of 
the two floes, which indicates a sticking probability less 
than 1 between floes. In particular, the sticking probabil- 
ity varies as different collision pairs since it is associated 
with the ratio of the eombinable region area to the total 
surface area of either floe in a pair. Therefore, interfa- 
cial tension-induced flocculation dynamics may result in 
a growth mechanism ruled by the relative sticking prob- 
ability during settling. 

It is worthwhile to note that the pressure p should be 
understood as an average effect resulting from different 
solid surfaces contributing to the external surface of the 
floe, which can be seen from Eq. (8). In other words, for 
different solid surface elements on the external surface of 
the floe, the pressure pi would be fluctuating around the 
average p deriving from Eq. (8). 



IV. THE MAXIMAL SIZE OF THE FLOC 

The interfacial tension-induced mechanism of floccula- 
tion implies that there is a maximal diameter of the floe, 
since size of the floe becoming large will abate the magni- 
tude of the additional pressure of the surface, which char- 
acterizes an ability of the floe's capturing free particles 
or floes. When sedimentation dominates over Brownian 
motion, the floe sweeps up smaller floes underneath and 
grows faster, and then, when can the floe's size attain the 
maximum? 

Let us assume that the floe is a sphere, then during set- 
tling, the upper hemisphere will receive few smaller floes 

due to gravity, conversely the lower hemisphere can com- 
bine large numbers of smaller floes. And we flnd that the 



pressure required to equilibrate a primary particle stuck 
on the floe, which is exerted by the interface between the 
external surface of the floe and the liquid, is minimum 
when the primary particle is enveloped on the bottom of 
the sphere (See Fig. 1). Assuming the actual additional 
pressure, thus the actual pressure acting on the primary 
particle, is averagely the same down the external surface 
of the floe and is determined by the size of the floe, we 
may know that the floe can not stick other floes together 
any more if it is not able to adhere to the primary parti- 
cle on its bottom. In other words, for the floe of certain 
a size, the nearer the particle or floe approaches the bot- 
tom of the main floe, the easier it is fixed, hence shape 
of the final floe generally is of cobble or strip. The case 
of the minimum pressure, corresponding to that of the 
largest floe, is exactly what we focus on. 

If enough particles are supplied to make the floe 
largest, then for a primary particle attached critically 
to the bottom of the floe, mechanical equilibrium gives: 
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where p' is the pressure acting on the primary particle 
from the interface, G is net weight of the primary particle 
in water. Therefore the maximal diameter Dm of the floe 
can be derived from Eq. (10) and Eq. (11): 
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Here we apply G = cb{ps — PL)gDQ = csApgD^, where 
CB is the shape factor, and As = ^tt-Dq. It is clearly seen 
that the maximal size Dm of the floe increases in a form 
of power law as the interfacial tension 7 grows, while it 

drops with the size Dq of the primary particle since the 
fractal dimension d/ is generally greater than 1 and less 
than 3. 

The expression of the maximal size Eq. (12) is a dy- 
namical result which excludes some geometrical factors 
such as competition between the floes in particles and 
spatial distribution of particles or floes. In other words, 
only if such dynamical parameters as 7, Do ete. can sat- 
isfy the critical condition of balance, the maximal size 
of the floe may be uniquely determined. However, ac- 
tual average size of floes in the system can not attain 
the maximal size because of clustering competition, spa- 
tial distribution, and sufficiency of particles. Therefore, 
the maximal size Dm is an imaginary quantity dependent 
only on dynamical factors. 



V. THE TENSION-DEPENDENT SETTLING 
VELOCITY 

In order to examine the effect of the interfacial tension 
on settling velocity of sediment floes, according to Eq. 
(1) and Eq. (7), we write: 



U=^^Dl- 



18 



J-/ 71 



(13) 



4 



the main floe 




the primary 
particle 



water 



FIG. 1: Schematic representation for mechanical equilibrium 
under the situation of the minimum additional pressure, thus 
the minimum p'. 



ing sediment system, since the pps indicates the essential 
factors affecting flocculation such as density of sediment 
hoc and binding force inducing flocculation etc., while 
electrolyte concentration is only a reflection of influences 
of solid adsorption action on flocculation. However, floc- 
culation factor expression Eq. (15) shows that F varies 
inversely as Ap and g. It is not surprising since smaller 
Apg will lead to larger average size of floes when bind- 
ing force is fixed, yet this still needs further experimental 
proof. 



where A = ^' d'l — = ^ df-i ■ 

We actually provide a general method for investigating 
settling velocity of sediment floes, since the dimensionless 
number A may remove the effect of the dynamical factors 
such as the interfacial tension and the primary particle 
size, and is associated only with geometrical factors like 
sediment concentration and spatial distribution of par- 
ticles, provided the growth of the average size {D'^f~^) 
indeed dynamically results from those reasons which lead 
to the the maximal size of the floe. In other words, the 
reason why the average size of floes cannot reach the max- 
imal could be that multi-body competition in particles 
under certain conditions reduces chances of the growth 
of each floe, which exactly arises from geometrical fac- 
tors. 

Then the average velocity of sediment floes can be ex- 
pressed through replacing Dm in Eq. (13) with Eq. (12) 
as follows: 
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We thus derive a setting velocity expression of sediment 
floes including dynamical information, which shows that 
the average settling velocity of sediment floes increases as 
solid-liquid interfacial tension becomes large in a form of 
power law, and has the decreasing relation with the pri- 
mary particle size, since Vladimir Nikora et al. 0| have 
proved that the fractal dimension df can affect relation- 
ship between velocity and size of the floe only when it is 
greater than 2. In addition, because of high electrolyte 
concentrations resulting in larger surface tensions, which 
is correct for inorganic salt solutes from the knowledge 
of physical chemistry, setting velocity would also increase 
with the electrolyte concentration. We may also obtain 
the expression of flocculation factor F: 
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which is deflned as F = U/Uq, where Uq denotes the 
classical Stokes relationship for solid spherical particles, 

namely Uo ~ Js^d'lT'' obvious that the remarkable 
effect on flocculation does not lie in electrolyte concen- 
tration but the size Dq of the primary particle compris- 



VI. COMPARISON WITH EXPERIMENTAL 
DATA 

A. The salinity-dependent settling velocity 

The effect of electrolyte concentrations on the settling 
behaviors of sediment floes may be derived from the rela- 
tionship between the surface tension and electrolyte con- 
centration C. The relationship between 7 and C is ex- 
hibited in the Gibbs adsorption equation: 



r2,i = - 
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where r2,i is surface excess on the Gibbs' dividing sur- 
face, C2 is solute equilibrium concentration, R is mole gas 
constant and T temperature. Therefore, if the isother- 
mal adsorption property r2,i(C2) of the interface is given, 
then the relationship between 7 and C2 may be obtained 
through integrating Eq. (16). The solid-liquid inter- 
face possesses more complex adsorption properties and 
accordingly experimental investigations on the relation- 
ship between F and C could be more reliable to study 
the salinity-dependent settling behavior of sediment floes. 
Nevertheless, many physicochemical experiments of sed- 
iment and soil show that the adsorption curve of some 
solid-liquid interfaces can be expressed by Langmuir ad- 
sorption isotherm, namely 



r2,i = 6- 



l + kC2 



(17) 



where b is the saturation adsorptive capacity, and k is 
adsorption constant. 

However, for the interface between solid and liquid, we 
are more interested in surface force relative to the solid 
rather than relative to the bulk liquid. It is known that 
the addition of salt to water will result in larger gas-liquid 
surface tension. This is because attraction of the salt ion 
in the bulk liquid for water molecules can draw more 
water molecules inside the bulk liquid so that more work 
would be needed to add a surface area, which shows that 
the adsorption quantity of the dividing surface relative 
to the bulk liquid is negative. In the same light, the 
adsorption of the solid to inorganic ions will also give rise 
to a negative adsorption quantity of the dividing surface 
relative to the bulk part between the dividing surface and 
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solid surface. Although the magnitude of the adsorption 
relative respectively to two different bulk parts, the bulk 
part between the dividing surface and solid surface, and 
the bulk liquid, is the same (the dividing surfaces in the 
two cases are different), the two quantities have opposite 
sign with each other, and this distinction will play a key 
role in determining properties of the interface tension. 

Thus, combing Eq. (16) with Eq. (17) and noticing 
sign of the adsorption quantity r2.i, we get 
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After the integration of Eq. (18), a tension-concentration 
relationship can be obtained: 



7 = fei?rin(l + fcC2) + 7o. 



(19) 



Here 70 denotes the interface tension between solid and 
pure solvent. Finally, the salinity-dependent settling ve- 
locity U and flocculation factor F can be expressed as 
follows, 
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respectively. It is apparent that both settling velocity 
of sediment floes and flocculation factor increase with 
increasing electrolyte concentration. However, real ad- 
sorption of solid in the solution is quite complex, even 
transitions between positive and negative adsorption of- 
ten occurs for some cases. Therefore it would be more re- 
liable that experimental methods determine the relation- 
ship between 7 and C to predict the effects of electrolyte 
concentration on settling velocity of sediment floes. The 
7-C form of Eq. (19) only reflects a simple case. In ad- 
dition, average size of floes dependent on the maximal 
size Dm can not be smaller than the pps Dqj hence the 
flocculation factor F is always greater than or equal to 
1, which indicates a cutoff concentration Cc which guar- 
antees {D) = XDjn > Dq. The cutoff concentration Cc 
actually is the minimum salt concentration causing floc- 
culation of particles of the Dq size in the present 7-C 
form. If 7 is a decreasing function of C, then Cc which 
satisfies F — 1 should represent a minimum salt concen- 
tration giving rise to maximal dispersion of floes. 




5 10 15 20 25 30 35 40 45 

The electrolyte concentration S (1/1000) 

FIG. 2: The fitting of experimental data to the model Eq. 
(23). The experiments were performed with coastal silt at the 
Bohai Bay in China of sediment concentration 4~5 kg/cm^, 
and the experimental temperature was about 293K. Solid lines 
denote fitting curves and Do** is the cumulative average size 
of the salt-free silt samples. 



B. The fitting of experimental data 

In order to examine the rationality of the assumption 
of tension-induced flocculation, we fit the typical exper- 
imental data which come from the report on sedimen- 
tation experiments of the coastal silt at the Bohai Bay 
in China by Wu Deyi to the electrolyte concentration- 
dependent sedimentation model Eq. (21) ll^,!!^- Wu 
Deyi et al. utilize the sedimentation balance to measure 
the electrolyte (A^a+) concentration-dependent settling 
velocity in still water of silt fiocs after analyzing the par- 



ticle size of salt-free silt sample by means of particle size 
analyzer. The salt-free silt sample can be gained through 
washing silt with distilled water and stirring and dispers- 
ing them in distilled water. They report an empirical 
formula between settling velocity and the salinity as well 
as the pps as follows, 

^ = h- Dq^-' -{S/lOOOf ^ (22) 

where Dq is the average size of the salt-free silt sample, 
S'/IOOO is the mass ratio concentration of Na~^, and the 
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TABLE I: Comparison of Fitted Values of F with Experimental Data as well as Values of Related Parameters 
of Best Fit 



Salinity F" 
[1/1000] (experimental) 


F 

(fitted) 


7-^ 


— 2/im 


%Difference 
between the two 
Fs 


B(best fit) 


fo(best fit) A;(best fit) 
[Mol/cm^] [L/Mol] 


1.3 
9.3 
12.9 
39.7 


32.4 
65.4 
113.5 
146.0 


26.74 
85.04 
97.69 
146.18 






17.5% 
23.1% 
13.9% 
0.12% 


1.79 X 10"^' 


0.99 62.49 


1.3 
9.3 
12.9 
39.7 


63.4 
129.0 
208.0 
258.0 


57.62 
159.3 
180.36 
259.89 


Do50 = 


1.35/im 


9.1% 
19.02% 
13.3% 
0.73% 


8.42 X 10"^^ 


0.997 98.8 


1.3 
9.3 
12.9 
39.7 


260.0 
400.0 
700.0 
600.0 


285.42 
488.93 
526.29 
661.50 


Dolo = 


0.55^im 


8.9% 
18.2% 
24.8% 
9.3% 


7.53 X 10~^^ 


0.998 1634.75 


ap — Us65 Us50 g^jjjj UslO 

ent Do**s. 


, corresponding respectively to difTer- 











numerical value of h is estimated as 4x10 ^. Due to 

F ^ Dq ^ , we get the fractal dimension df in an 
average sense df ~ 2.2 which is empirically reasonable 
according to a large number of results of experimental 
investigations on the fractal dimension of dispersions. 
Let us rewrite Eq. (21) as 

F = B- [bRTHl + fc^) +7o]^(Ap5)^i^7^, 

(23) 

and take R = 8.31 x lO^(dyne-cm-mol-iK-i), T = 
293(K), 70 = 72.92(dyne-cm-i), Ap = 1.5(g-cm-3), and 
g = 1000(cm-s^^). Here the unit of k is [L/mol] and B is 
a dimensionless parameter. We use the least-squares fit 
of the model Eq. (23) to experimental data reflecting the 
relationship between F and S to estimate values of the 
parameters B, &, and k. The results are shown in Fig. 2 
and Table. I. 

It can be seen from the results of fit that the floc- 
culation factors obtained from the model Eq. (23) are 
close to those from real experiments. Therefore the as- 
sumption of the tension-induced flocculation of sediment 
particles are reasonable. The more experimental data 
would attain better fits, nevertheless, percentage differ- 
ence between Fs obtained by fits and experiments may 
still provide a verification of the rationality of the model 
proposed. In particular, the magnitude of the fitted pa- 



rameter k which possesses definite physical sense accords 
with our understanding of the adsorption of metal ions 
on sediment particles which is generated from lar ge n um- 
bers of experiments, for example, Ref [2^ l30l l3lL 132 . l33| . 
etc. In addition, the fact that values of the fitted pa- 
rameter b in different measures are approximately equal 
indicates that the parameter b may indeed represent such 
a physical quantity as saturation adsorptive capacity. 



VII. CONCLUSION 

This paper proposes a theoretic model which attributes 
dynamical reason of flocculation of floes to the interfacial 
tension between solid particles and liquid, and provides 
a general method of analyzing sedimentation of sediment 
fiocs. The resulting expression of settling velocity of sedi- 
ment floes can give a good description of the relationship 
between flocculation factor and the electrolyte concentra- 
tion without loss of physical sense of fitted parameters. 
In fact, actual interaction is so complex that more me- 
chanical factors need to consider. For example, during 
aggregation ruled by Brownian motion while interacting 
fiocs being effectively small, other forces have to be intro- 
duced into balance conditions. Further study should be 
aimed at direct confirmation of the present model with 
real experiments. 
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